Replica theory for fluctuations of the activation barriers in glassy systems 



OO ■ 

o ■ 
o : 

(N . 

Oh: 

m ; 

m ■ 
(N ; 

^\ 
? ■ 

C/3 \ 



I 

o 



> 

00 
00 
On 

m 

ON 

o 

OO 

o 



Maxim Dzero/ Jorg Schmalian,^ and Peter G. Wolynes'^ 

Center j or Materials Theory, Rutgers University, Piscataway, NJ 08854, USA 
^Department of Physics and Astronomy and Ames Laboratory, Iowa State University, Ames, I A 50011 
'^Department of Chemistry and Biochemistry and Department of Physics, 
University of California, San Diego, La Jolla, CA 92093 
(Dated: September 23, 2008) 

We consider tlie problem of slow activation dynamics in glassy systems undergoing a random first 
order phase transition. Using an effective potential approach to supercooled liquids, we determine 
the spectrum of activation barriers for entropic droplets. We demonstrate that fluctuations of the 
configurational entropy and of the liquid glass surface tension are crucial to achieve an understanding 
of the barrier fluctuations in glassy systems and thus are ultimatively responsible for the broad 
spectrum of excitations and heterogeneous dynamics in glasses. In particular we derive a relation 
between the length scale for dynamic heterogeneity and the related barrier fluctuations. Diluted 
entropic droplets are shown to have a Gaussian distribution of barriers, strongly suggesting that 
non-Gaussian behavior results from droplet-droplet interactions. 



I. INTRODUCTION 

The glass transition and glassy behavior are dynami- 
cal phenomena characterized by slow relaxations and a 
broad spectrum of excitations. Recently strong exper- 
imental evidence for dynamical heterogeneity in glassy 
systems, i.e. for spatially varying characteristic time 
scales for relaxation, has been obtained through NMR^i 
and nanometer-scale probing of dielectric fluctuationsi. 
The wide distribution of excitations, as seen in bulk di- 
electric measuments, is therefore believed to be caused 
largely by spatial variation of the characteristic time 
scales; see also Ref."*. Regions in a glassy system sep- 
arated by only a few nanometers relax on time scales 
different by many orders of magnitude. 

Whether these purely dynamical phenomena can be 
explained in terms of the underlying energy landscape 
is still debated. The question has been positively an- 
swered on the level of the mean field theory of glasses. 
The dynamical, ideal mode coupling theory^ and energy 
landscape based replica mean field theories were demon- 
strated to describe the same underlying physics, yet from 
rather different perspectives^'''. The idealized version of 
the mode coupling theory"" falls out of equilibrium below 
a temperature Ta > Tg, where Tg is the laboratory glass 
temperature (see below). This is consistent with static, 
energy landscape based mean field theories that find 
an emergence of exponentially many metastable states, 
Afms exp {const. V), at the same temperature Ta- Here 
V is the volume of the system. Thus, the configura- 
tional entropy Sc = ks \ogMms becomes extensive and 
only vanishes at the Kauzmann tempeature Tk < Ta- 
Unfortunately, neither the emergence of slow activated 
dynamics, nor the possibility of spatial heterogeneity are 
correctly captured within mean field theory, making it 
impossible to use mean field theory alone to directly ex- 
plain the heterogeneous dynamics observed in numerous 
glass forming liquids. On the other hand, mode coupling 
theory at higher temperatures (T > Ta) has been shown 
to give a reasonably accurate account of the dynamics of 



liquids. At lower temperatures {T < Ta), bulk thermo- 
dynamic properties of glassforming liquids, most notably 
the temperature dependence of the configurational en- 
tropy and the related Kauzmann parados^, agree with 
the results of the energy landscape based mean field the- 
ory (see Refs.— liS). This motivated the development of 
a theory for glassforming liquids that has its foundation 
in the mean field theory, but that does take into account 
activated dynamics and spatial heterrogeneity which are 
beyond the strict mean field approacb^iS.. The resulting 
random first order (RFOT) theory of glasses demon- 
strates that ergodicity above the Kauzmann temperature 
Tk can be restored via "entropic droplets" , in which a 
region of the liquid is replaced by any of an exponentially 
large number of alternatives. The entropy of the possible 
alternatives acts as a driving force for structure change. 
This driving force is off-set by the free energy cost of 
matching two alternative structures at their boundaries. 
This conflict gives a free energy barrier for activated 
motions^i^. Various consequences of entropic droplets 
have been analyzed in Refs.— li^i^ii^ using for concrete 
calculation a density functional description of the liquid 
state. These quantitative results have been found to be 
in good agreement with experiment. Within RFOT, the 
number of metastable states determines the activation 
barrier of the glassy state. While this is similar to the 
Adam-Gibbs theory for glassesi^, there exist important 
distinctions between the two approaches, see Ref.^-'^. In 
agreement with the RFOT approach but not with the 
assumptions of Adams-Gibbs argument the complexity 
of a correlated region increases as the glass transition 
is approachedi^. The need for instanton like events to 
understand the long time dynamics in glasses is also sup- 
ported by recent findings extended mode-coupling theo- 
ries for dense fluids^*. 

More recently, an effective Landau theory for glasses, 
based on the replica method of Ref.-'^^, has been devel- 
oped. This framework naturally allows for activated dy- 
namics beyond mean field theory2£i2i. This approach 
based on instantons offers a formal justification of the 
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entropic droplet approach using the rephca approach to 
glasses. It reproduces and thus confirms several results 
of the random first order transition (RFOT) theory of 
structural glasses^. In both approaches, a mean droplet 
activation energy Fi occurs that is determined by the 
configurational entropy density Sc (for details see below) : 



Ft 



(1) 



p2i- vanishes linearly at the Kauzmann temper- 



Sc OC 

ature Tk < Tg. The exponent v relates the droplet size 
R and the configurational entropy density: R oc s~'^. d 
is the spatial dimension. Below we assume that d — S. 
In its most elementary version the replica Landau the- 
ory yields ly = 1— i^. A renormalization of the droplet 
interface due to wetting of intermediate states on the 
droplet surface^ was shown to yield = 2/d, leading to 
Fi cx Ts~^ and correspondingly to a Vogel-Fulcher law 
^ ^^P3/^ ^ for the mean relaxation time 



To exp I 



Toexp 



Ft 



(2) 



A softening of the surface tension due to replica symme- 
try breaking of the instanton solution was also obtained 
in Ref;^^ and seems to be a first correction containing the 
effects that lead to a reduction of the exponent from 
the value v — I. 

Numerous experiments on supercooled liquids not only 
depend on the mean barrier, but are sensitive to 
the entire, broad excitation spectrum in glasses^. Most 
notably, the broad peaks in the imaginary part of the 
dielectric function e" (uj) are most natually understood 
in terms of a distribution g (r) of relaxation times, such 
that 



(3) 



e (w) cx / drgir) -3. 

1 + {ujt) 



Similarly dynamical heterogeneity with spatially fluctu- 
ating relaxation times yields nonexponential (frequently 
streched exponential) relaxation of the correlation func- 
tion 



<t>{t)= drgir) e-"\ 



(4) 



and allows for an explanation of nonresonant hole burn- 
ing experiments^^, even though interesting interpreta- 
tions of the non-exponentiality based on " dynamical ho- 
mogeneity" exist as well2^. Other effects that are most 
likely caused by a distribution of relaxation rates include 
the breakdown of the Stokes-Einstein relation D 



between the diffusion coefficient Z? of a particle of size L 
and the viscosity r f^i'^^ . 

These experiments all call for a more detailed analysis 
of the fluctuations 



of the activation barriers and, more generally, of the dis- 
tribution function p [Ft") of barriers. The latter yields 
the distribution function of the relaxation times 



9ir)^p{Ft)'^ 



(5) 



through r [F^j — tq exp (^-^f^ ■ For example, in case of 

a Gaussian distribution of barriers one obtains a broad, 
log-normal distribution of relaxation rates: 



5(r) 



: exp 



log^ (t/t) 
2A 



with 



x = SFty{kBTY 



(6) 



(7) 



and T from Eql2l While the distribution, Eq[6l does not 
yield precisely a stretched exponential form for the cor- 
relation function, it can often be approximated by 



</.(t)c.exp {-(t/rf) 



(8) 



with (3 = (1 + X) '. Furthermore, the study of higher 
order moments of p (^Ft^ is important to determine 
whether the distribution is indeed Gaussian or more com- 
plicated. 

In this paper we determine the fluctuations SFt^ as 
well as higher moments of the distribution of the barrier 
distribution function p {Fty We start from the replica 
Landau theory of Ref and use the replica instanton 
theory of Refsi ^°'^^ . We then generalize the replica for- 
malism to determine barrier fluctuations in addition to 
the determination of the mean barrier Ft . Our approach 
is an important step to gain insight into the whole spec- 
trum of activated processes in glassy systems in general. 

The role of barrier fluctuations within the RFOT the- 
ory of structural glasses was first discussed by Xia and 
Wolynes in Ref.i^.. Variations in the configurational en- 
tropy density were assumed to be the main cause of 
the fiuctuations of the barrier. Fig. [TJ Starting from 
Ft OC sl~'^'^, barrier and configurational entropy density 
fiuctuations are related by 6Ft cx s~'^''6sc, which yields 
SFf^ (X sJ^'^'^Ss'^. Standard fluctuation theory then de- 
termines the entropy fluctuations in a region of size R'^ 
as Ssl (X AcpR~'^ where the configurational heat capac- 
ity remains finite as T —^ Tk (for a detailed discussion of 
Acp see below) . This reasoning finally yields 



5Ff^ 



OC s, 



-du 



OC R"^. 



(9) 



In the approach of Ref.— , barrier fluctuations diverge 
with the volume of the entropic droplet. Our results will 
demonstrate that the leading droplet size dependence of 
8Ft'^ agrees with Eql9] and Refi^. We also show that 
there are additional terms originating in the fluctuations 
of the interface tension of entropic droplets that con- 
tribute to bFf^ . In addition we analyze the temperature 
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FIG. 1: Appearance of an exponentially large number of 
meta-stable configurations and the possibility for a system to 
realize these configurations gives rise to the configurational 
entropy. The latter serves as the driving force for the struc- 
ture change in the glassy phase. The energy cost to match the 
boundaries between the different density configurations gives 
rise to the energy barriers. The fluctuations of the configura- 
tional entropy give rise to the energy barrier fluctuations. 



dependence of the static barrier fluctuations and of the 
exponent [3 for stretched exponential relaxation as well 
as the variation of the dielectric response with frequency 
and temperature. 

This paper is organized as follows. In the next sec- 
tion we introduce the replica effective potential formal- 
ism to analyze barrier fluctuations in glassforming liq- 
uids. We motivate the approach by starting from a den- 
sity functional approach to liquids and using an equi- 
librium replica theory to derive the mean field theory as 
starting point of our calculation. Next we analyze fluctu- 
ations of the conflgurational entropy in bulk, summarize 
the effective potential approach of Ref?l£ and the instan- 
ton calculation that yields the mean barrier F^. Finally 
we demonstrate how higher moments of the barrier dis- 
tribution can be derived within the same formalism. At 
the end of the section we present our results for physical 
observables. The paper is concluded with a summarizing 
section. 
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FIG. 2: The typical energy landscape in the structural glassy 
phase as a function of density (left panel) and the effective 
potential, which is used to describe the transition between 
the nonergodic glassy states and the ergodic liquid state (right 
panel). 



dom hard sphere packing22., instead of the periodic crys- 
tal lattice positions. The free energy of this amorphous 
solid was shown to have a global minimum at a = and 
a local minimum for finite a. If a ~ y-^/a _^ ^j^g 
particles are delocalized and the system in an ergodic 
liquid state with homogeneous density po = ^. Finite 
a corresponds to a frozen amorphous solid i.e. a glassy 
state. For T > Tk the amorphous solid is metastable 
with respect to the liquid, and higher in free energy by 
TSc- It is of course always metastable with respect to the 
crystalline solid. Essentially the same mean fleld physics 
can be described in terms of a formally more precise 
replica approach, introduced by Monasson"^*^. Here one 
determines the partition function in the presence of a bias 
conflguration p (r): 



Z[p]= / Dpe 



,-l3<P[p]-g ! d-'x(p(x)-p{x)f 



(10) 



Here, j Dp... corresponds to the statistical sum over all 
density configurations of the system. The free energy of 
a bias conflguration is 



II. THE REPLICA EFFECTIVE POTENTIAL 
FORMALISM 

A. Motivation and cloned replica approach 

We start our description of glassy systems from the 
point of view that there exists a density functional, (p [p], 
that describes a supercooled liquid undergoing a mean 
field glass transition. Initially, such an approach was 
used in Ref.-^ where it was shown that it allows one to 
describe the emergence of a metastable amorphous solid. 
Fig. [21 Following the classical approach to freezing into 
ordered crystalline statea^^, the density was assumed to 
be p(v) = (f )'^^^ Here, a determines the 

Lindemann length, a~^/^ over which particles are local- 
ized. In distinction to crystallization, the mean positions, 
Ti, of an amorphous solid were taken to be those of a ran- 



f[p\ = -T log Z[p\. (11) 

Physically / \p\ can be interpreted as the free energy for 
a metastable amorphous configuration of atoms, for ex- 
ample with density p{x) being a sum over Gaussians dis- 
cussed above. In the replica formalism, no specific config- 
uration like this needs to be specified in order to perform 
the calculation. Rather, the assumption is made that the 
probability distribution for metastable configurations is 
determined by / [p\ according 

P[p]cxexp(-/3eff/[p]) (12) 

and is characterized by the effective temperature Toff = 
P^i > T. This allows one to determine the mean free 
energy 

F = y Dpf[p]P[p] (13) 
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and the corresponding mean configurational entropy 



Sc = -J DpP[p]\ogP[p\ 



(14) 



The approach of Ref.- - was successfuhy used to de- 
velop a mean field theory for glass formation in super- 
cooled liquids^ yielding results in detailed agreement 
with earlier, non-replica approaches^!. As shown by 
Monasson^, the mean values F and Sc can be deter- 
mined from a replicated partition function 

Z [m] — J £)™pe~''^"=i '^^'''''+^^""''=1 

_ ^ (15) 
via F ^ -§^^mF (m) and S'c = (m) with 



T 

F (to) = log Z (m) 



(16) 



and replica index m = . It follows from the replicated 
free energy F(Toff): 



F = 

s: = 



d {F (Teff ) /Teff ) 



dF{T,s) 
dT,s ' 



(17) 



These results are in analogy to the usual thermodynamic 
relations between free energy (F — > F(Tcs)), internal 
energy {U — > F), entropy {S Sc) and temperature 
(T Tcs), see also Ref<^. If the liquid gets frozen 
in one of the many metastable states, the system cannot 
anymore realize its configurational entropy, i.e. the mean 
free energy of frozen states is F, higher by TSc if com- 
pared to the equilibrium free energy of the liquid. Fig. 

m 

If the replica solution is taken to be marginally stable, 
so that the lowest eigenvalue of the fluctuation spectrum 
beyond mean field solution vanishes, it has been shown^i 
that Teff agrees with the result obtained from the gen- 
eralized fluctuation-dissipation theorem in the dynamic 
description of mean field glasses^. Typically, the as- 
sumption of marginality is appropriate for early times 
right after a rapid quench from high temperatures. In 
this case Teff > T for T < Ta, i.e. the distribution func- 
tion of the metastable states is not in equilibrium on the 
time scales where mode coupling theory or the require- 
ment for marginal stability applies. Above the Kauz- 
mann temperature it is however possible also to consider 
the solution where T^g = T, i.e. where the distriburion 
of metastable states has equilibrated with the external 
heat bath of the system. Since we are interested in a 
system where ergodicity is restored for Tk < T < Ta 
we use Teff = T. Technically this is a supercooled liquid 
and not a non-equilibrium aging glass. Below the Kauz- 
mann temperature this assumption cannot be made any 
longer (at least within the mean field theory) as it leads 
to a negative configurational entropy, inconsistent with 
the definition EgfTH 



The physically intuitive analogy between effective tem- 
perature, mean energy, F and configurational entropy to 
thermodynamic relations, Eq[T7l suggest that one should 
analyze the corresponding configurational heat capacity 



Cn — Tci 



dm 



(18) 



Using the distribution function Eq ll2l we find, as ex- 
pected, that Cc is a measure of the fluctuations of the 
energy and configurational entropy of glassy states. We 
obtain for the configurational heat capacity 



{spy 



{sScY 



(19) 



where {SFy ^ F^ ~F and {6Sc) ^ S^-Sc • Here the 
mean values F and Sc are determined by Eq |17l The fiuc- 
tuations of the configurational entropy and frozen state 
energy are then determined by S*^ = J DpP [p\ log^ P [p] 

and F'^ = J DpP [p] f [p] , respectively. Both quantities 
can be expressed within the replica formalism in terms 
of a second derivative of F (to) with respect to to. For 
example it follows: 



dr 



;mF (m) = -^ (f^ - F^ 



(20) 



It is then easy to show that EqlTH] holds. With the intro- 
duction of the configurational heat capacity into the for- 
malism we have a measure for the fluctuation of the num- 
ber of available metastable states from its mean value. 
The analogy of these results to the usual fluctuation the- 
ory of thermodynamic variables^ further suggests that 
Cc also determines fluctuations of the effective tempera- 
ture with mean square deviation: 



(21) 



Since Cc is extensive, fluctuations of intensive variables, 
like Toff, or densities, like Sc = Sc/V, vanish for infinite 
systems. However, they become relevant if one considers 
finite subsystems or small droplets. In the context of 
glasses this aspect was first discussed by Donth^^-. 

Using standard techniques of many body theory we 
can proceed by expressing the replica free energy F (m) 
as a functional of the liquid correlation functions of the 
problem (see e.g. Ref.— ): 

H[X] = ^ (Tr (xo'x) + TTr In X + $ [x]) ■ (22) 



Here the mean field result of F (to) is determined by the 
stationary point of the effective Hamiltonian: ^ = 
with replicated density-density correlation function 



Xab [X, x') = {Spa (x) dpb {x')) 



(23) 



The Luttinger-Ward functional [x\ is determined by the 
nature of the density-density interaction and determines 
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the self energy S = — 12 , see also Refsi^i^. For example, 
expanding the density functional, </> (p) ~ 4'{pQ + 5p), 
into a Taylor series in 5p, yields to lowest order in 5p: 



X^Jx,x') (24) 



where v-}, = ^ rfps""* ■ what follows we will formulate 
the theory in terms of the collective variable Xah (x, x'), 
instead of the original density fluctuations and write 

Xah (x, x') = X (X, x') bah + Cab (x, x') (25) 

where Cab = if a = 6. Cafc(x, x') depends on two 
spatial variables x and x' . We simplify the problem by 

assuming that its Fourier transformation, Cab { ^^^ ; i 

with respect to the relative coordinates x — x' can be 
written as 



Cab (r; q) = qab (r) paS (q) . 



(26) 



with liquid structure factor, (q) and r = (x + x')/2 
is the center of mass coordinate. The important col- 
lective variable of our theory is qab (r), which plays the 
role of a spatially varying Debye- Waller factor. Within 
mean field theory we expect below the temperature Ta 
that qab = 9* (1 — 5ab) is replica symmetric with q* of 
order unity, while q* vanishes in the equilibrium liquid 
state for T > Ta- Physically the Debye- Waller factor q* 
contains the same information as does the localization 
parameter discussed above. Finally, it is important 
to stress that a replica symmetric approach of qab in the 
present formalism is equivalent to one step replica sym- 
metry breaking in the conventional replica language^*^. 

In order to keep our calculation transparent we will not 
determine H [q] from an explicit microscopic calculation 
for supercooled liquids. It was demonstrated by Franz 
and Parisi^ that this is possible. However, in order be 
able to make simple calculations beyond the mean field 
theory explicitly including droplet fluctuations, we start 
from a simpler Landau theory in the same universality 
class^^^: 



H 



with 



bqbcqca 



(27) 



h[qab]^^{^qabY 



t 2 
2lab 



-q!^, + lqL (28) 



and replica index a, 6 = 1, . . . , m. ao is a length scale of 
the order of the first peak in the radial distribution func- 
tion of the liquid and Eg is a typical energy of the problem 
that determines the absolute value of Ts^. In addition 
the problem is determined by the dimensionless variables 
t, u, w and y, which are in principle all temperature de- 
pendent. We assume that the primary T-dependence is 



that of the quadratic term, where t = ^-^^ . Making this 
assumption, along with the polynomial form of the La- 
grangian, requires constraints between the other Landau 
parameters to ensure thermodynamic consistency from 
the fluctuation formulas. In Appendix A we give esti- 
mates for the parameters of Eqsl27| and [28] obtained 
from fits to experimental data for o-terphenyl (OTP), 
a well studied fragile glass forming material. Formally, 
Eqsl27] and [28] can be motivated as the Taylor expansion 
of Eq[22] together with the functional $ of Eq[M] In 
what follows we further simplify the notation and mea- 
sure all energies in units of Eq and all length scales in 
units of qq. 

The mean field analysis of this model Hamiltonian is 
straightforward. Inserting a replica symmetric ansatz 
qab = 9* (1 ~ Sab) into H and minimizing w.r.t q* yields 
q* = or 



2y 



(29) 



Nontrivial solutions exits for t < tA = with q* (<a) = 
which determines the mode coupling temperature 



Ta = tAEo + To. Inserting q* of Eql 
for the replica free energy: 



I into H [q] yields 



F{m) 
V{m,-1) 



t 



,*2 



w -\- u (m — 1) 



yq*\ 



(30) 



The mean configurational entropy density, as determined 
by ^ = Tv^l^ ("^)L^i is gi^en by 



T~c - ^q*' - ^q 



*3 



(31) 



Inserting q* of Eq[2H] yields the result that vanishes at 

9y 



tj^ = ^ with q* {tx) = j^- Close to Ik it follows that 



tx ( t — tK\ T — Tk 

oc 



4y V Tk 



T 



K 



(32) 



as expected. At tA one finds TaSc (Ta) — jUp-- 

We next determine the configurational heat capacity 
of the mean field theory discussed above. From Eq|5Dlwe 
obtain 



^ _ 2Vm ft^^^ 



w-M(2-3m) 3 y^^^4 
Toff V 2" 3 ^4 



(33) 



If Teff = T (corresponding again to the equilibrium be- 
havior between T^ and Ta) it follows that Cc (Tk) = 

^1 (ly ) ' s° iii^SLt comparing it with the expression for 
the configurational entropy (1311) at t ~ we can write 



(34) 



Notice that w/Autfc must equal 1/Tfc for complete ther- 
modynamic consistency. For temperatures close to but 
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FIG. 3: Effective potential fl[pc] as a function of the overlap 
at the fixed temperature Tk < T < Ta- At the Kauzmann 
temperature an entropy crisis takes place with vanishing en- 
tropic advantage of the liquid compared to the frozen solid. 



above Tk the energy fluctuations decrease as 

8u + w t — 
2 ~ 



Ce ~ Cc {Ik) 1 - 



(35) 



As the temperature approaches Ta the configurational 
heat capacity decreases faster than linearly. It reaches a 
value Cc {Ia) = ^ ^"gg^ ^ signalling that the system be- 
comes unstable close to Ia for w > 8u. This is analogous 
to the higher order replica symmetry breaking suggested 
by Tarzia and Moored on the basis of the same Landau 
action. Cc {t — > Ia) approaches Cc {Ia) from above with 
an infinite slope: 



Cc {t ^ Ia) 

Cc (tA) 



24 



1/2 



8u 



(36) 



The full temperature dependences of the mean field 
overlap q* and configurational heat capacity are shown 
in Fig. [3] While q* or did not depend on the parame- 
ter u of the Hamiltonian, Eqsl2l] and[28l the configura- 
tional heat capacity does. We may therefore consider u 
as a parameter that determines the strength of entropy 
fluctuations. 



study the behavior for arbitrary values of the overlap 
between two states, including inhomogeneous droplet so- 
lutions. This is possible within the effective potential ap- 
proach introduced in Refsji^i^i^. We use this technique 
to formulate our replica instanton and barrier fluctuation 
theory. 

We are interested in the regime above Tk (i.e. for 
t > Ik) and consider the partition function of a system 
with constrained overlap between the particle density 
configurations pi (r) and p2 (r): 



n^(Pc(r)-p2(r)pi(r)). (37) 



The configuration p2 is assumed to be in equilibrium 
and unaffected by pi. Averaging the free energy, 
—Tlog Zp^,p2, weighted with the equilibrium probability 
for p2 yields the corresponding mean energy for the over- 
lap Pc-. 



fl [pc 



-T 



(38) 



Here Z is the equilibrium partition function Z = 
J Dpe~^*[''l . Using a replica approach SI [pj can be writ- 
ten asi^ 

n [p^] = ~T lim lim — f - l) , (39) 

ri^O m^O TO V / 



with 

m 

X n S{pc{r)-piir)pt{r)). (40) 

r,b=l 

il [pc] can be determined by analyzing an n + to times 
replicated problem with additional constraint for the 
overlap between certain replicas. In complete analogy to 
the previous paragraph one can introduce a field theory 
of the overlap, but now with order parameter 



r p 
T 

p q 



(41) 



B. The effective potential approach 

In order to investigate the transition state for the de- 
cay of a specific metastable frozen state, it is necessary to 
have a technique that agrees with the mean field replica 
approach discussed in the previous section in case of a 
homogeneous solution q = q*, but that allows one to 



where r is an n x n matrix, g an m x m matrix and the 
nxm matrix p must obey the additional constrained that 
Pib (r) = Pc (r) for 'ib=i,...,m- We have the same effective 
Hamiltonian, Eqsl27l and [28l onlv with q replaced by Q, 
i.e. 

/m 
DQe-''"^Q^Y[Sipi,-pc). (42) 
6=1 
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Following Franz and ParisilS we use rap 
Pab = ^a.iPc and it follows 

Hp. M = [ d'^ri 2mh \pc]+Y.h [qab 

\ ab 



and and 



oZ^ lab 



MbcQca - Upi X! 
abc be 



(43) 



The remaining g^ft-depcndent problem is formaly similar 
to the original one of EqsHT] and [551 but in an external 
field upl- 

Using these results we find that the effective potential 
can be written as 



with 



n [pc 



dr 



(44) 



(45) 



We start with a homogeneous replica symmetric ansatz 
for Qab = q (Sab — 1) with equal off diagonal elements, 
q, and zero diagonal elements. It was demonstrated 
that this replica symmetric ansatz is again equivalent to 
one step replica symmetry breaking in the conventional 
replica language^. We take the limit m ^ and find 



%M = (^q) + 2h [p,) + uplq \uq\ 
Vm 3 



(46) 



If we extremize this with respect to q, i.e. perform the in- 
tegration with respect to q at the level of a homogeneous 
saddle point approximation, it follows that 



tq — {w — u) q^ + yq^ 



UPc 



(47) 



Depending on the value for the overlap pc, q is in general 
different from the value q* obtained within the replica 
approach of the previous section. If we however insert 
q back into Hp^ (q) and determine the effective potential 



n[pc] 

for Pc 



- and for p, 



Pc (9 (Pc)) we find that it has minima 







q{p^) = q\ with q* of EqlH 
in complete agreement with the approach of the previ- 
ous section. The value of the effective potential for the 
metastable minimum pc = q* is identical to the mean 
configurational entropy discussed above 



(48) 



In order to analyze the local stability of the solution at 
the minimum of f2 (pc) we determine the replicon eigen- 
value of the problem. Following Ref."*^ gives 



A,. 



t — 2wq + 3yq 



(49) 



for the lowest eigenvalue. Inserting q from Eq. II47P yields 
close to Ik a positive eigenvalue: ~ i/f — 5 (i — i/f). 



0.75 
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FIG. 4: Mean field value for the configurational lieat capacity 
(main plot) and overlap (inset) for y — 1.82, ui — 2.72, u = 
0.385 (see text) as a function of reduced temperature {t — 
tK)/{tA — tK)- Exponentially many states emerge below Ia 
with the discontinuous jump of the overlap. 



while the replicon eigenvalue vanishes at Ia as A^ — 
2tA ( tj^ ] • This behavior is completely consistent 

with the results obtained in Ref.''* for the random Potts 
model. The homogeneous solution is marginal at Ta and 
stable below Ta- Except for temperatures close to T^, 
the mean field approach is locally stable with respect 
to small fiuctuations. The effective potential approach 
is not restricted to overlaps Pc that minimize 17 [pc] ■ In 
general q {pc) is given by the solution of Eg 1471 Inserting 
this solution into Hp^ (q) gives (pc)- The result is shown 
as broken line in Fig.®. It turns out that for interme- 
diate values of Pc, between pc = and Pc = q* an ad- 
ditional replica symmetry breaking of qab occurs^. For 
the present model this effect was also studied in Refi^. 
This is related to the well known effect that an external 
field can cause replica symmetry breaking above Tx^, 
yielding one step replica symmetry breaking for qab- The 
result for fl {pc) shown in Fig. 2] has been obtained, in- 
cluding this additional replica symmetry breaking. Phys- 
ically this additional replica symmetry breaking as Pc 
" goes over the hill" was discussed in Ref^Si. We also note 
that the generic behavior for fl (pc) as shown in Fig. [His 
very similar to Ts^(g*) of Eq[3Il if we simply plot this 
function for arbitrary q* . As required both expressions 
exactly agree at saddle points and minima. 

In summary, we showed that on the level of spatially 
homogeneous mean field solutions, the effective potential 
approach of RefsJ^i^i^ yields results in complete agree- 
ment with the physically transparent approach of Refi^. 
Since it allows for arbitrary overlap between two states, 
we use it as starting point for the determination of acti- 
vated events within an instanton calculation. 
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C. Instanton theory for activated events 

At the mean field level a glass at T < Ta is frozen in 
metastable states and characterized by the overlap q* of 
Eg between configurations at distant times. For tem- 
peratures T > Tk above the Kauzmann temperature, 
the free energy of such a frozen state is higher by the 
configurational entropy TSc compared to the liquid state 
that is characterized by q* = 0. Thus, for Tk <T <Ta 
the mean field glass is locally stable ( the replicon eigen- 
value Xr is positive) but globally unstable with respect to 
the ergodic liquid. The effective potential il (pc) shown in 
Fig. 3 suggests that the decay modes for the frozen state 
are droplet excitations, similar to the nucleation of an 
unstable phase close to a first order transition. This sit- 
uation was analyzed in Refj ^°'^^ . Using this approach, 
an analysis of the effective droplet size for a Kac-type 
model was performed in Ref4^ In agreement with the 
RFOT theory^, the driving force for nucleation is the 
configurational entropy, leading to the notion of entropic 
droplets. 

Instanton solutions for entropic droplets are deter- 
mined from 



10 



5pc (r) 







(50) 



where we allow for spatial variations of the overlap Pc (r). 
Recently, Franz has given a transparent dynamical in- 
terpretation for the solutions of EqlSOP. Assuming 
that the effective potential is characterized by spatially 
anisotropic but replica symmetric solutions qab (r) = 
q (r) {5ab — 1), we perform the integral over the qab in 
Eql45l at the saddle point level. We find a solution 
q (r) = pc (r) with 



V^p, (r) 



dVjPc (r)) 
dpc (r) 



where 



V{Pc)^\pl~^pl + \pt 



(51) 



(52) 



EalSlI admits an exact solution in the thin wall limit 



th( ^ 



- th ( - 



(53) 

where the integration constant zq is a function of i, w and 
y. R is the droplet radius and ^ is the interface width 
given by 



4an 



y/3y{2q* - q*^Y - 6tK + At 



(54) 



The final expression for the function Z[){t) is more in- 
volved and is given in Appendix C. 

Inserting the solution Eg 1531 into the expression for the 
effective potential EqUJ] we calculate the value of the 
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FIG. 5: We show the temperature dependence of the ratio ^/_R 
[R is the size of the droplet) together with the temperature 
dependence of interface width ^(t) (inset). 



mean barrier. The latter is determined by optimizing 
the energy gain due to creation of a droplet and energy 
loss due to the surface formation. As a result for the 
mean barrier we find (reintroducing the energy scale Eq 
and length scale aq) 



32^ao 2 



The droplet radius 



R 



Mai 



(55) 



(56) 



is determined from the balance between the interface 
tension and the entropic driving force for nucleation. 
The barrier energy of Eq l55l determines the relax- 
ation time T of EqlSl Furthermore, q^ = q*{t = Ik) 
is the order parameter at the Kauzmann temperature. 
When temperature approaches the Tk the radius of the 
droplet as well as the mean barrier diverge. One finds 
lim Ft (X [t — tK)^'^ and lim R oc (t — tK)^^- Since the 

t^tK t—ftK 

droplet interface ^ remains finite as i — > the thin wall 
approximation is well justified close to the Kauzmann 
temperature. On the other hand, R and ^ become com- 
parable for temperatures close to Ta and the thin wall 
approximation breaks down. We see that replica Landau 
functional calculation predicts a rather diffuse droplet 
near the laboratory Tg. Combining i? cx {t ~ tx)^^ and 

cx {t — tpc), we obtain v = 1 for the exponent that re- 
lates the droplet size R and the configurational entropy 
density: R a s~'^ . 

In Ref . [s^ it was demonstrated that for a similar model 
the interface free energy is exponentially small for a large 
system and in any finite dimension the one step replica 
symmetry breaking state does not exist. As was already 
pointed out in Ref. [s^ , our approach is very different in 
scope and in its conclusions. While Ref. is concerned 
with the absence of replica symmetry breaking in the 



ultimate equilibrium state, our conclusions are relevant 
for the nonequilibrium sitation. Since we obtain a finite 
barrier height of the frozen mean field solution, our result 
offers a mechanism for equilibration on time scales < 3> t 
demonstrating that our approach and the conclusions of 
Ref. m are consistent. 

D. barrier fluctuations 

To analyze the barrier fluctuations we start from the 
barrier for a given density configuration p : 

P'pcP ~ -^PcP ^ -^q' ,p (57) 

where Fp^ p — —T log Zp^^p is determined by the con- 
strained partition function of Eg 1371 and Fq, ,p the corre- 
sponding energy with homogeneous overlap, i.e. Fq'^p = 
—Tlog Zp^j^j,^oo)—tq' .p- When we analyze variations of 
activation barriers we need to keep in mind that both, 
ground state and transition state contributions to F| „ 

^ Pc tP 

are statistically fluctuating and are correlated in general. 
Fluctuations in the first term of Eq |57l correspond to vari- 
ations of the free energy of the localized instanton. On 
the other hand, fluctuations in the second term corre- 
spond to variations of the homogeneous background. In 
principle cancellations between both occur that are prop- 
erly included in the analysis of F^^ ^ defined above. 

Barrier fluctuations are then characterized by the sec- 
ond moment 



SFt^ = F^ip - Fl^p (58) 

where the average is, just as for the analysis of the mean 
barrier, with respect to the density configuration p. The 
second term in Eql58lis the square of the mean activation 
barrier F'^ and was determined in the previous section. 
Thus, we can concentrate on first term. Using Eg 1571 
it follows that the first term in EqlSS] consists of three 
terms: 



Fpt,p - F^^ p + Fq, p - 2Fp^^pFq*^p (59) 

In what follows we consider these three contributions sep- 
arately. The first two terms correspond to independent 
fluctuations of the droplet and the homogeneous back- 
ground, while the last term measures their mutual corre- 
lations. We find that this last term in Eql59] is propor- 
tional to the surface area of the instanton. Correlations 
between droplet and homogeneous background terms re- 
sult from their mutual interface. 

The detailed analysis of the three contributions to Fp'^^p 
of Eg 1591 is summarized in appendix B. In what follows 



we summarize the results of this derivation. For the first 
two terms in EqinHit follows that 



F^',P = Fq,J-2T I d'rh[q*] 

Fi; = Fp::^ ~2T j <fvh[p, (r)]. (60) 

The first expression gives fluctuations of the configura- 
tional entropy for the homogeneous problem, while the 
second one describes the energy fluctuations of the con- 
figurations with spatially heterogeneous overlap. 

The result for the homogeneous problem, F^, p can al- 
ternatively be obtained from our analysis of the configu- 
rational heat capacity given in Eg |19l Using the replica 
formulation of Eql^D] together with the explicit result of 
EqEOl yields 

F^~f'^ ^2VT - + fg^^) (61) 

where V is a volume. If we now recall the definition of 
h[q] Eql2Hl one readily observes that the first equation in 
Eg 1^ coincides with Eg. (pT]) . 

The derivation of the third contribution Fp^ pFqt p to 

i^p^,p in Eqinni which describes the energy correlations 
between the instanton and environment, is also per- 
formed in the Appendix B. The result is: 

Fp,,pFq,,p = F^^p X F^Tp -T J {2h[pc (r)] 

+upl{r) - uq*p,{r)) . (62) 

Combining Eqs (101 with EqlB^l we finally obtain for the 
second moment of the barrier fluctuations: 



^-^^ jd\ (MpJ - h[q*] + upliv) - uq*pl{v)) . 

(63) 

This result for the second moment of the barrier fluctua- 
tions depends on the value of the homogeneous overlap q* 
of Eql29]as well as the instanton solution Pc(x). The thin 
wall limit result for Pc(x) is given in Eq[531 Inserting 
these expressions into the Eg l63l vields 

= A{R^ + pR") , (64) 

where R is the radius of the droplet. In the thin- wall 
approximation, the explicit expressions for the coefficient 
A is 
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where the value of qi is computed using an effective po- 
tential V{qi) = V{q*), i.e. it is the turning point of the 
intstanton, see Fig. 3. Deriving ([65|) we also made the 
following choice for the parameters q*{tK) = iw/iy = 1 
(see Appendix A for more details). This choice is of par- 
ticular convenience as it allows one to express the result- 
ing expressions in terms of the ratio t/tj^ and u. The 
length scale p{t) can be compactly as follows: 



Pit) 



A 



2uq* 



C2(t) 



Au 



C3(i) +C4(t) 



(66) 



where Cn{t) are functions of the ratio t/tK only and are 
defined via 



c„(i)= / [g™-p^(z)]d^ 



(67) 




(t-t^)/(t^-t^) 



FIG. 6: Temperature dependence of the length scale p{t) de- 
scribing the contribution to the barrier fluctuations due to 
surface energy fluctuations of the structural droplets. We 
also show the temperature dependence of the ratio p/R {R is 
the size of the droplet). 



We have computed these integrals numerically for n = 
2, 3 and 4. Temperature dependence of the surface length 
scale p{t) is shown on Fig. [HI 

Before we proceed with the analysis of our results, it 
is instructive to write down the expressions for the co- 
efficient A and p{t) for temperatures close to Ik- When 
t ~ Ik, one readily finds qi ~ txit — tK)/y, so that ex- 
panding expression in the square brackets (j65p up to the 
linear order in t — we have A{t) ~ {kBTEQ/aQ)[AQ + 
Ai-it- tK)] + 0{{t ~ tKf) with 



87r 
= — 



Ai = 47r ( — 



1 



(68) 



where once again we employed q*{tK) = 1 (see Appendix 
A) . Thus we find that close to Ik the bulk contribution to 
the barrier fluctuations is governed solely by the value of 
the phenomenological parameter u. Expression for p[t ~ 
tx) = PK + PiK'it — tK) can be derived straightforwardly 
using our results above together with the expansion for 
the functions c„(i) and £,{t) ~ + 5{t — tK)/t%'^ ■ 

Here we provide the expression for the value of p{t = tx)'- 



pit 



9t 



K 



1/2 
K 



2u 

CiK - C3K + — I C2K - -C3K 



(69) 



where we have used the notations Cnx = CnitK)- We 
have found that the coefficients CnK have the following 
values: C2k — 1.25, C3/f ~ 1.18 and c^k — 1.11. 

Thus, in addition to the bulk term, which behaves in 
a similar way as the formulation in Refj^^ (see EqlHl) , the 
replica theory yields a surface term, resulting from fluc- 
tuations in the interface. In the inset of FiglS] we show 
the T-dependence of the length scale pit), demonstrating 
that the surface term becomes gradually more important 
as the temperature increases. The value of pitx) remains 
non-singlar, so that the volume term is dominant. This is 
consistent with the view that the droplet interface close 




(t-t^)/(t^-t^) 



FIG. 7: Temperature dependence of the energy barrier fluc- 
tuations, scale pit) (bottom inset). The latter describes the 
contribution to the barrier fluctuations due to surface energy 
fluctuations of the structural droplets. 



to Tk is smaller then the droplet radius, while both are 
comparable as T approaches Ta, similar to the behavior 
close to a spinodal'*^. The resulting temperature depen- 
dence of p{t) is shown on Fig. [HI 6F^'^ is shown in Fig. [7] 
and the exponent of the strcchcd exponential relaxation 
(3c shown in Fig. [S] 

The plots are constructed by setting the parameters 
of the theory to re-produce experimentally relevant val- 
ues of Tx, Tg and as well as configurational entropy 
above the glass transition using o— terphenyl as the ex- 
ample of the glass transition (see Appendix A for details). 
There is one free parameter left to be fixed, either y or 
w = 3?//2. We have determined the value of y in the com- 
parison with the experimental value of the mean barrier 
at T = Tg given by 7^ = TgsdTg) with the theoretically 
derived expression ((SS]) . This procedure yields y — 1.82, 
w = 2.73 and u = 0.385. As we can see from Fig. [51 
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(t-t^)/(t^-t^) 



FIG. 8: Temperature dependence of the exponent /So gov- 
erned by the fluctuations of the isolated droplets together 
with exponent /3 corrected by the effects of the droplet-droplet 
interactions. 



PciTg) — 0.12 , a value significantly reduced compared 
to the exponential behavior where f3 — 1. This value 
for Pq is lower than the experimentally found values of 
PexpiT — Tg) ~ 0.5. This discrepancy may signal that 
the present approach overestimates the strength of the 
energy barrier fluctuations an effect that may be related 
with non-Gaussian fluctuations of the activation barriers. 
It is however clear that these may be described as instan- 
ton interaction effects as discussed by Xia and Wolynea^^. 
Our analysis yields a Gaussian distribution of barriers 



Pg (pi) = 



1 



exp 



/ 



V 



(^Fi-Fty 

2SFt^ 



(70) 



which yields the following estimatefor the exponent of 
the streched exponential relaxation 



f3c 



i + sFty{kBTy 



-1/2 



(71) 



Our results for the temperature dependence of Pg are 
shown in Fig[5] The calculation that led to this result was 
based on the assumption of single instanton events, i.e. 
entropic droplets were assumed to be diluted. Obviously 
this is not addressing the fact that distict droplets in- 
teract. The mosaic picture underlying the RFOT theory 
is clearly based upon the view that such droplet-droplet 
interactions occur and are in fact crucial. The impact of 
droplet-droplet interactions for the distribution of barri- 
ers was analyzed in Ref^^^. Here it was pointed out that, 
as soon as the droplet size becomes larger than the size 
i? of a mean droplet (Eq.55 in our theory), boundary ef- 
fects will limit the size of a droplet, leading to a cut off 
of the distribution function. This leads to the modified 
distribution function 



P{Ft) 



5 (f* - Ftj 



pi < Ft 
Ft >Fi 



(72) 



which leads to a reduction of the mean square deviation 
of the barriers by a factor 1/4, compared to the Gaussian 
distribution. This corrects the exponent of the streched 
exponential relaxation to 



1- 



i + -SFty{kBTy 



-1/2 



(73) 



where SFi"^ is still the result, Eq.63, of the dilute droplet 
calculation. In FigElwe compare Pg and p. Obviously, 
the reduction of the barrier fluctuation width due to 
droplet-droplet coupling leads to an increase of the ex- 
ponent p. We mention however the resulting value of P 
at the glass transition is P{Tg) ~ 0.23. This suggests 
the present Landau functional indeed over-estimates the 
fluctuation effects. This probably implies that the dif- 
fuseness of the droplet from the Landau functional is also 
overestimated. The fluctuations do vanish at Ta, consis- 
tent with the emergence of simple exponential dynamics 
at a significantly high temperature. We did not attempt 
to vary the parameters of our model to reach perfect 
agreement with experiments for OTP, as the simplified 
Landau functional is clearly an oversimplification of the 
physics of real glasses. It does however demonstrate the 
generic and qualitative trends for barrier fluctuations in 
the RFOT theory. 

Finally, in Fig. [3] we show the dielectric finction 
e" {uj, T) as finction of frequency for different tempera- 
tures. The peak in e"{uj,T) dramatically broadens as 
temperature approaches the glass transition tempera- 
ture, Tg, as predicted from the static barrier distribu- 
tions. For this figure we used the second moment SF^^ of 




500 
CO/O) (T) 

win ' 



2.5e+05 



FIG. 9: Frequency dependence of the imaginary part of the 
dielectric susceptibility above the glass transition. The broad- 
ening of the peak at the glass transition signals an appearance 
of wide spectrum of excitations. 

the barrier distribution function and then assumed that 
p {Pt^ is Gaussian. In case of a Gaussian distribution. 
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all moments can be expressed in terms of the second mo- 
ment 

SFik ^ (jFi^y^^ (74) 

where /^^ = for k odd and fi^ = 2'T (^) /\/tt for 
k even. In particular it holds for higher order moments 
= 0, = 3, fif = 0, and = 15 etc. 
We can explicitly analyze higher moments of this dis- 
tribution and thereby demonstrate explicitly that p (-f"'") 
is Gaussian. To be precise, we have only shown this 
for the first six moments but strongly suspect this to be 
true in general. The calculation of higher moments are 
quite cumbersome but they formally are a straightfor- 
ward generalization of the method we used to determine 
the second moment SFi'^. In case of the k-th moment 
the block structure of the matrix Q of Eq |B4l consists 
of (fc + 1) X (fc + 1) blocks where the total dimension of 

Q is (^n + J2i=i "^fe) X (^n + "^fe) • The rather te- 

dious matrix algebra is then most easily analyzed using 
computer algebra software, such as Mathematica. Thus, 
we only list the results here. We analyzed all moments 

= J (^F^ -F^y p{F^) dF^ (75) 

up to fc = 6 of . We find 

SFt^ = 0. (76) 
In case of the fourth moment we analyze: 



6Fi^ = iFt[p]-Ft)^ 

= Fti [p2] - 4Ft^ [p2] Ft+3Ft + 6Ft- 7^(77) 

it follows 

Wt* = 3[2V [pc-, qo] - [Pc] - nho^{q*)f (78) 

where Tihom and Timh are defined below Eq. IB8I Evalu- 
ating this expression yields 

jF^^?,{6F^y (79) 

We further find that the fifth moment vanishes as well 

= (80) 
and that the sixth moment is given as 

Jf^ ^ lb {p^y . (81) 

Thus, up to the sixth moment the barrier distribution 
of droplets is Gaussian. Thus, dilute entropic droplets 
should indeed have a Gaussian distribution of barriers, 
strongly suggesting that the observed non-Gaussian be- 
havior of the effective barrier distribution results from 
droplet-droplet interactions. 



III. SUMMARY 

In summary, we have shown how the energy barrier 
fiuctuations SF^"^ as well as higher moments of the static 
free energy barrier distribution function p {F^^ can be 
computed using the replica Landau theory of Ref.JS and 
use the replica instanton theory of Refsi^i^i . We have 
generalized the replica formalism to determine barrier 
fluctuations in addition to the det ermin ation of the mean 
barrier Fi. Barrier fluctuations 5F'^'^ consist of a dom- 
inating bulk contribution at temperatures close to the 
Kauzmann temperature along with a contribution stem- 
ming from fluctuations in the interface. The latter dom- 
inate the barrier fiuctuations well above the Kauzmann 
temperature. Fluctuations of the surface energy enter 
prominently in the recent work of Biroli et alj^ on di- 
rectly computing point-to-set correlations in liquids. It 
is very nice to see that the effect they found naturally 
emerges from the replica instanton calculation frame- 
work. It is interesting to speculate that these surface 
energy fluctuations could account for the observed sys- 
tem specific deviations from the Xia-Wolynes prediction 
of a direct relation between the stretching exponent and 
ACt,. Chemical trends in these deviations may provide 
clues on this. JFt^ increases as the typical droplet size 
increases, i.e. as the temperature is lowered towards Tr-. 
Clearly, our theory only applies until the system falls out 
of equilibrium at the laboratory glass temperature Tg 
with Tg > Tk- In the framework of the present instan- 
ton approach, the energy barrier distribution function 
turns out to be Gaussian. The skewness of the observed 
relaxation time distribution is, in our view, an effect due 
to the interaction of spatially overlapping instantons. It 
can be captured, at least in spirit, by an extended mode 
coupling approach, incorporating instantons. 
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APPENDIX A: ESTIMATE OF THE MODEL 
PARAMETERS FOR O-TERPHENYL 

The Landau functional is quite convenient for for- 
mal use but there are variety of ways that it may be 
mapped onto real fiuids. In this appendix we show how 
the model parameters of the Landau replica potential. 
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Eqs[27| and [2ll can be found by fitting the results of 
the mean field theory to experimental data. We use 
the concrete example o-terphenyl (OTP), a well known 
glassforming material. Values for the effective dynamical 
freezing temperature and the Kauzmann temperature of 
OTP are Ta = 285X^4 and Tk = 202.7K.^2i59. On the 
other hand, a typical value for the glass temperature isS& 
Tg ~ 243...246K. Here, the glass transition temperature 
is the temperature at which the viscosity reaches a value 
of 10"'^^g/(cm s), which is also where the mean relaxation 
time reaches values of lO^s. The melting temperature of 
OTP is r,„ = 329K with an entropy jump at melting of 
AS*™ = 52.28^^^. 

The configurational entropy at Tg and Ta are Sc {Tg) — 



21.5- 



2. 59R and {Ta) = 39- 



'molK ~ ^. "^-i^ c^ii^i ^c\--^AI — •-"^inolK ~ 4. 69R, 

respecticely^i^, where R is the molar gas constant. To 
obtain the entropy per spherical object or bead one has 
to divide these results by ub — 3.7, the number of 
beads as defined by the procedure in Ref.^^. This yields 
Sc {Tg) = = 0.7R as well ws {Ta) = 1. 26i?. 



The heat capacity change at the glass transition ii 



.55 



Acp {Tg) = 111.27— J— - 0.48^ (Al) 
molK gK 

where we used Ig = 4.31 x 10~^mol. The specific and 
molar volume of the system are^?: 

cm^ cm'^ 

v{Ta) = 0.918 = 213.2 

g mol 



cm 



cm 



v{T„) = 0.893 = 207.3 -. 

g mol 



(A2) 



The typical volume per particle V/N — is then char- 
acterized by the length lo- It follows V/N = 344. 23A^, 
such that lo = A.35A. In comparison, the van der Waals 
radius of OTP was given as^ ''vdw — 5.7 A. 

If we introduce the energy and length scales into the 
problem, the mean field theory results obtained from the 
Hamiltonian EqsHT] and are 



qo {Tk) 
Ta 
Tk 

TaSc {Ta) 
Sc (T ~ Tk) 
TkCc {Tk) 



2w 

V Tq Ay 



To 1 



3 ° 192y3 
47r^3 Eqw / 2w\^ t — tfc 

47r gfi'o 2u f2wY 



3 %[J 3 



(A3) 



The order parameter at T^ should be large of order unity. 
We chose qo {Tk) = 1 which gives w = |?/. Using the 



relation Ta — Tk = -^^Eo this yields 

Eqw = 1975. 2K = 9. 74Tk 



(A4) 



If we further make the reasonable choice of oq = 
0.87rydw for the length scale oq, Sc {Ta) of Eq lA3l gives 
precisely the experimental value listed above. We can 
also determine the change in heat capacity Acp {Tg) from 
the same parameters. Eq. I A4I yields the result that 



t-tK 

tK 



0.31- 



T~T 



K 



Tk 



(A5) 



so that we obtain for the temperature dependence of the 
configurational entropy: 



Sc (T) ~ Acp {Tg) 



T-Tk 
T„ 



(A6) 



where the mean field result for the slope is Acp {Tg) = 
4.36i?. This compares reasonably well with the experi- 
menal value Acp (Tg) = = 3. 62R. Note, how- 

ever, that Acp as given by Eq |A6l is different from the 
configurational heat capacity Cc, Eq. ((33| . It holds 
Sc ~ ^Cc{tK) i.e. f^Cc{tK)0.253 = Acp (T,,), 

i.e. both are comparable provided u ~ O.lui, a value 
we will assume in what follows thus ensuring thermo- 
dynamic consistency from the two fluctuation formulas. 
In summary, the existing thermodynamic measurements 
can determine several parameters and allow a consistent 
description of several thermodynamic data. Still, one free 
parameter, either w oi y remains since only the ratio w to 
y is determined. We use this freedom to obtain the cor- 
rect value of the mean barrier at the glass temperature. 
Finally, we note that the above result also imply that 
= 0.06 6 i.e. about half of ^-^=^ = 0.125. It is 
worth mentioning that the values of the Landau parame- 
ters are consistent with the random first order framework 
as assumed. Finally, we summarize our results for the 



TABLE I: Parameters for the Landau functional 



^dvW 



0.1 1.125 1.066 2.73 0.385 1.82 



model parameters by fitting the experimentally relevant 
quantities to the data obtained for OTP in Table [J 



APPENDIX B: DETAILS OF THE REPLICA 
CALCULATION FOR BARRIER 
FLUCTUATIONS 

In this appendix we summarize the derivation of the 
various contributions to the barrier fluctuation given in 
Eg 1551 We start our calculation with the third term in 
Eg 1591 The other two terms are simpler to determine 
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and can be obtained as specific limits of the third one. It 
holds: 



^Pc,p2 ^q* ,P2 



T- 



JVpe-P^^P^ logZp^,plogZ,.,p 



lim (3^12-3^1-3^2) 

™l,2— >0 



(Bl) 



where we introduced 

3^12 = T 

yi = 



2 J Va,P q 



Z 



m\miZ 



3^2 = 



m2 



mim2 /Ppe-^^W 



(B2) 



where Z = J Vpe /^^W. To analyze the first term in 
Ea lBll we use again the replica trick and write: 



X12 = lim 3^12 

■mi 2—^0 



lim 



I?(^l 



— lim lim ^"''^ 

mi ,7712 — ^0 n — >0 



-/3H[(^i]^ml ^m2 , ^ 

mim2 (/I?(^e-/3«M)""^ 

(n;mi ,7712) 



17111712 



(B3) 



Here we introduced a quantity Xp"'^.'^'™^'' that can be 
determined from an n + mi+m2 replicated problem with 
order parameter Q: 



/ 



Q 



r 



p q 

u ip 



(B4) 



It holds 



-7-(n;mi,m2) 
Pc;9o 



PQe-^«[Qin'5(Pc(x)-pi.(x)) 

x,a 

l[S {qo ~ qU^)) , (B5) 



i/[Q] has the same form as Easl27landl28lif one uses the 
matrix Q of Eg 1B4I instead of the original replica variable 
q. The sub-matrices of Q are given as follows: r is an 
n X n matrix, p is an n x nii matrix, qo is an n x m2 
matrix, u is an nii x mi matrix, v is an m2 x m2 matrix, 
and ijj is an nii x m2 matrix. The matrices p and qo obey 
the additional constraints: 



pib(x) = pc(x) for V = 1, ...,mi; 
qia = go for V a = 1, ...,m2. 



(B6) 



which is enforced through the (5-function above. We as- 
sume for a replica symmetric instanton solution that and 
the elements of V'a/S are all the same: 



V'a/3(x) = ■0(x). 



(B7) 



Similary we assume Vaf3 = (1 — Sai3)v. 

The analysis of H[Q,] tedious but straightforward and 
yields that it can be written as a sum of three terms: 

H[Q] miHinh[pc;u] + m2Tihom.[qo;v] 

+mim2V[V']. (B8) 

Here 

nho^[q*;v] = j d'^{2h[q*] + {m2-l)h[v] 

- (to2 - 1) vq*'^ - i (rn2 - 1) 
refers to the homogeneous problem. In addition, 

'H^nh[Pc■,u] = J d^x{2h[pc] + (m^i - l)h[u] 

-(mi - l)u(x)p^(x) 



--(mi - l)(mi - 2)u3(x) 



is the contribution of the inhomogeneous instanton. Fi- 
nally, the coupling term between the two is given as 

Vlij] = J d^x{2h[i;] + (mi - l)uo(x) 

+ (m2 - l)w]V^2(x) - 2p,(x)g*V(x)}. 

In order to determine the various matrix elements of Q 
defined in Eq |B4l we perform saddle point approxima- 
tions with respect to the numerous variables. We find 
u(x) = r(x) which has been calculated earlier. Similarly 
it follows V = q* . The saddle point equation for ip is then 

Pc{^)qo - -VVXx)-FtV'(x)-aV'^x)-hyV^(x()B9) 
+ [ro(x) +q*]^(x) 

land it follows that the solution of this equation is 
'0(x) = Pc(x). With these results we are in the position 
to calculate the barrier fiuctuations. 

First we re- write expression Eq.(Bl) by employing the 
saddle-point solution as follows: 



lim 



1 



1 -e" 



-)n2/3Hftom[9*,9*] 



mi, 2^0 mim2 
_g-mi/3Hi„fc[pc;ro.i'l]_|_ 

^-mil3'Hinh[Pa\ro,ri\-m2l3'Hhom[q' -mim2/3V[pc] ^ (BIO) 

Expanding the exponents up to the second order in m^'s 
and taking ro.i = Pc yields the expression (57) in the 
text. Equations (55) can be obtained using the same 
steps which lead us to (jBlOp . For example, we have: 



F2 

Pc,p2 



lim ( 

11,2-^0 mim2 V 



1 - e 



-m2/3'Hi„h[pc;i"o,'"i] 



-mi/3Wi„h[pc;ro,ri] 



-{rni+m2)l3'Hinh[Pc\'ro,rx]-mirn2pVi„h[pS^ (BH) 
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where we introduced Vinh[Pc] = J (P^h[pc]- Upon taking 
the hmit mi 2 — > and using the rephca symmetry rg = 
ri = Pc we recover the second expression in Eqi59l 



In terms of these functions the expression for zq (t) reads: 



APPENDIX C: EXPRESSION FOR THE 
INTEGRATION CONSTANT zo 



zo{t) = - log 



(C2) 



In this Section we provide an expression for the inte- 
gration constant zq which enters into Eq. (|53p for the 
interface profile. First, we define the following functions: 



qoit) 



2y 



1 - 



8t 



K 



t — - t, 



3qlit) - 1 



(CI) 



Finally, we remark that the integral, which determines 
the length scale p{t) (|66|) can be computed exactly using 
relations above. The resulting expression is quite cum- 
bersome and will not be listed here. 
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